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Unconventional  Energy  Sources  for  Ice  Control 
at  Lock  and  Dam  Installations 

TATSUAKl  NAKATO,  ROBERT  ETTEMA  AND  KEIICHI TODA 


INTRODUCTION 

A  major  concern  for  operating  a  navigation  lock  and 
dam  during  frigid  weather  is  ice  formation  along  lock 
walls,  at  lock  miter  gates  and  at  dam  gates.  Ice  collars 
create  numerous  operational  problems,  including  re¬ 
ducing  the  effective  lock  width,  creating  difficulty  in 
operating  miter  gates,  interfering  with  vessel  passage 
and  making  it  necessary  to  reduce  tow  width.  Gates 
made  inoperable  by  being  frozen  to  ice  sheets,  or  by  be¬ 
coming  encrusted  in  ice,  pose  potentially  adverse  con¬ 
sequences  for  flow  regulation.  It  is  therefore  extremely 
important  that  appropriate  methods  be  found  for  con¬ 
trolling  ice  growth  at  lock  and  dam  installations. 

Two  general  methods  can  be  adopted  for  controlling 
ice  growth:  mechanical  and  thermal.  Mechanical  meth¬ 
ods,  which  involve  removing  ice  that  has  already  formed, 
are  part  of  routine  operations  generally  followed  at 
locks  and  dams.  A  variety  of  devices  are  used  to  remove 
ice,  including  backhoe  scrapers,  chainsaw  cutters,  pike 
poles,  high-pressure  water  jets  and  steam  jets.  These 
methods,  however,  can  be  arduous  and  occasionally 
hazardous.  Thermal  methods  entail  two  alternatives:  a) 
applying  sufficient  waste  heat  available  from,  for  exam¬ 
ple,  a  powet  plant  to  keep  the  temperature  of  water  in  the 
lock  chamber  above  its  freezing  temperature;  or  b) 
keeping  lock  walls,  gates  and  pier  walls  adequately 
warm,  in  the  vicinity  where  ice  would  form,  such  that 
ice  growth  is  inhibited. 

Conventional  energy  sources,  primarily  electricity 
from  power-generation  utilities  or  hydrocarbon  fuels, 
have  proven  to  be  comparatively  expensive  when  used 
for  ice  control  at  lock  and  dam  installations.  Therefore, 
unconventional  energy  sources,  such  as  solar  energy, 
wind  energy  and  sensible  heat  from  groundwater,  ap¬ 
pear  attractive.  This  investigation  assesses  the  feasibil¬ 
ities  of  various  unconventional  energy  sources  for  ice 
control.  It  includes  a  review,  through  several  site  visits, 
of  ice-control  practices  generally  used  at  lock  and  dam 


installations.  Also  discussed  are  the  heating  require¬ 
ments  for  ice  control  at  locks  and  the  characteristics  of 
heat  transfer  through  concrete  lock  walls  or  pier  walls. 


ICE-GROWTH  PROBLEMS  AT 
NAVIGATION  LOCKS  AND  DAMS 

The  two  ice-growth  processes  of  concern  are  ice- 
collar  formation  along  lock  walls  and  lock  miter  gates, 
and  ice  growth  against  dam  gates.  Under  most  circum¬ 
stances,  ice-collar  growth  is  a  form  of  border  (or  shore- 
fast)  ice  growth,  forming  at  the  water-surface  elevation 
predominantly  maintained  in  a  lock  (Fig.  1 ).  This  may 
be  at  either  the  upstream  or  the  downstream  pool  eleva¬ 
tion,  depending  on  individual  lock  operation  practice. 
Or,  as  it  is  fairly  common  to  allow  some  continuous 
throughflow  of  water  to  inhibit  ice-cover  formation 
over  the  entire  lock  chamber,  an  ice  collar  may  form  at 
an  elevation  between  these  pool  elevations. 

Relatively  quiescent  water  adjoining  the  walls  of  a 
lock  experiences  less  vertical  mixing  of  flow  and  more 
rapid  cooling  of  the  flow  surface  than  the  water  toward 
the  center  of  the  lock.  Consequently,  as  is  the  case  with 
many  water  bodies  exposed  to  frigid  air,  ice  growth 
commences  at  fixed  boundaries  and  then  thickens  and 
progresses  toward  the  center  of  the  lock.  Raising  and 
lowering  the  water  level  within  the  lock,  together  with 
accumulation  of  snow,  frazil  and  ice  debris  and  abrasion 
during  barge  impact,  cause  ice  collars  to  become  some¬ 
what  rounded  in  cross  section.  Although  the  dimensions 
of  ice  collars  vary  with  growth  conditions,  typical  maxi¬ 
mum  dimensions  reported  by  Corps  of  Engineers  (COE) 
personnel  are  up  to  1  m  in  height  and  1  m  in  depth  (Fig. 
2).  Under  some  circumstances,  it  may  be  possible  for  ice 
growth  to  extend  down  much  of  the  lock  lift  height. 
Conceivably,  this  might  occur  if  a  series  of  collars  are 
allowed  to  form  because  the  water  surface  is  held  for 
significant  periods  at  several  elevations.  Alternatively, 
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Figure  I.  Ice-collar  growth  along 
lock  walls  of  Starved  Rock  Lock 
and  Dam,  Illinois  River,  January 
1982. 


Figure  2. Typical  geometry  and  maximum  dimensions  of  ice 
1  collars. 
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Figure  3 .  Icing  over  a  lock  wall 
below  an  ice  collar  at  Starved 
Rock  Lock  and  Dam,  Illinois 
River,  4  February  1978. 
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Figure  -4.  Ice  sheet  im¬ 
mediately  upstream 
from  L&D  No.  5  on 
the  Mississippi  River, 
23  January  1 986 . 


Figure  5.  Ta inter  gate 
w  ith  splash  and  spray 
i'  ing  at  L&D  No.  5  on 
the  Mississippi  River, 
26  January  1986. 


repeated  wetting  of  lock  walls  may  cause  ice  to  accrete 
as  an  icing  (Fig.  3).  Such  icing  of  lock  w  alls  is  usually 
less  problematical  than  ice-collar  grow  th. 

Ice-growth  problems  on  dam  gates  originate  from 
several  causes: 

•  Upstream  gate  faces  may  become  frozen  and  stuck 
to  adjoining  river  ice  covers  (Fig.  4); 

•  Ice  growth  due  to  w  ater  splash  and  spray ,  on  the  rear 
of  gates  as  well  as  areas  surrounding  tainter-gate 
trunnions  and  roller-gate  ring-gear  and  rack  devic¬ 
es.  may  render  them  inoperable  (Fig.  5). 


•  Gates  may  become  unbalanced  by  uneven  ice 
growth  such  that  they  are  difficult  to  operate:  and 

•  Gate  side  seals  are  prone  to  ice  accretion  because  of 
water  seepage  as  well  as  water  splash  and  spray. 

Sometimes,  due  to  the  effects  of  water-level  fluctua¬ 
tions.  a  crack  separates  the  ice  cover  from  ice  formed  as 
a  collar  along  the  gate  face  and  gate  piers  (Fig.  5). 

Ice  problems  at  locks  and  dams  and  in  navigable 
waterways  are  discussed  by  Zufelt  and  Calkins  ( 19X5 ). 
They  provided  a  summary  of  problem  severity  for  the 
Allegheny,  Illinois,  Mississippi  and  Ohio  risers. 


FIELD  PRACTICES  FOR  ICE  CONTROL 

To  review  current  ice-control  practices,  we  visited 
several  lock  and  dam  (L&D)  installations  during  Janu¬ 
ary  and  March  1986.  Lockmasters  were  interviewed 
regarding  ice-control  problems  and  practices.  Because 
of  the  limited  number  of  site  visits,  the  following  com¬ 
ments,  which  are  based  on  the  daily  experience  of  the 
lockmasters  interviewed,  are  not  necessarily  conclu¬ 
sive. 

We  visited  two  locks  and  dams  on  the  Mississippi 
River  (L&D  No.  5  at  Minnesota  City,  Minnesota,  and 
L&D  No.  6  at  Trempealeau.  Wisconsin).  Most  lock 
facilities  on  the  Mississippi  River  are  generally  closed 
when  thick  ice  forms  in  the  river.  However,  dam  gate 
structures  must  be  operational,  even  during  severe  win- 
termonths,  to  control  river  stages.  Therefore,  much  em¬ 
phasis  seemed  to  be  focused  on  how  to  keep  large  rol  ler 
gates  ice-free.  At  L&D  No.  5  a  large  carriage-mounted 
steamer  is  used  in  conjunction  with  steam  lances  to  de¬ 
ice  dam  gates.  Both  installations  use  electric  side-seal 
heaters  embedded  in  the  roller-gate  pier  walls.  Three 
side-seal  heating  elements  (2.3  kW  per  element)  that 
operate  together  are  placed  in  each  pier  wall  pocket 
(Fig.  6).  The  lockmasterof  L&D  No.  6said  that  de-icing 
with  side-seal  heaters  can  be  accomplished  within  about 
45  minutes  of  their  application.  Quartz  heaters  (2.5-kW 
units)  are  also  used  extensively  to  melt  ice  that  forms  on 
gate  racks.  Figure  7  shows  two  quartz  heaters  being  sus¬ 
pended  over  a  roller-gate  rack. 

The  power  requirements  for  controlling  ice  growth 
on  dam  gates  at  these  L&Ds  can  be  estimated  on  the 
basis  of  these  heater  capacities.  Consider,  for  example, 
L&D  No.  5,  which  is  equipped  with  six  6. 1-m-high, 
1 8.3-m-long  roller  gates.  If  two  quartz  heaters  (2  x  2.5 
=  5  kW)  and  three  side-seal  heaters  (3  x  2.3  =  6.9  kW) 
are  used  for  each  pier  wall,  a  total  of  1 42.8  k  W  would  be 
needed.  Because  three  or  four  roller  gates  are  generally 
sufficient  to  regulate  low  winter  river  flows,  only  half  of 
this  power  would  be  required  to  keep  three  roller  gates 
ice- free. 

At  Marseilles  L&D  (Marseilles,  Illinois)  on  the 
Illinois  River,  electrical  heating  cables  were  installed 
on  the  top  of  the  miter  gates  to  keep  the  gates  ice-free. 
Insulation  pads  were  placed  over  the  heating  area. 
Quartz  heaters  similar  to  those  used  on  the  Mississippi 
River  installations  were  tried  around  miter-gate  hinges. 
However,  they  were  found  to  be  ineffective.  Instead, 
high-pressure  water  jets  were  used  to  wash  solid  ice  off 
gate  quoin  areas.  The  reason  quartz  heaters  failed  to 
work  here  is  that  Marseilles  Lock  is  operational  during 
most  of  the  winter.  Therefore,  areas  surrounding  gate 
hinges  are  subject  to  frequent  changes  in  water-surface 


elevation,  which  would  produce  thicker  ice  than  in  the 
roller-gate  case,  where  ice  forms  more  slowly  from  a 
combination  of  spray  and  seepage  water.  Large-scale 
transverse  bubbler  systems  were  also  used  to  deflect  ice 
away  from  miter-gate  recesses.  The  bubbler  systems 
must  be  kept  running  continuously  to  avoid  pipe-freez¬ 
ing  problems.  Properly  designed  bubbler  systems  gen¬ 
erally  do  not  require  continuous  operation. 

Three  L&D  facilities  were  visited  on  the  Ohio  River. 
Emsworth  L&D.  located  at  Emsworth,  Pennsylvania, 
uses  bubbler  systems  extensively  in  miter-gate  recesses 
as  well  as  near  the  approach  and  exit  sills  of  the  lock 
chamber.  A  hydropower  generator  rated  at  300  kW  was 
newly  in  service  to  provide  the  power  for  general  opera¬ 
tion  and  for  ice  control.  The  utility  bill  without  turbine 
output  during  the  winter  was  about  four  times  larger 
than  that  with  the  operating  turbine.  Dashields  L&D  at 
Glenwillard.  Pennsylvania,  has  an  ungated  fixed-crest 
dam.  so  it  does  not  experience  gate-related  ice  prob¬ 
lems.  A  large  air-bubbler  system  with  a  0. 1 7-m-Vs  com¬ 
pressor  is  used  to  flush  ice  away  from  miter-gate  recess¬ 
es.  Kerosene-bumer-type  heaters,  with  capacities  of  44 
kWh  per  unit,  are  lowered  along  lock  walls  to  melt  ice 
collars  formed  there.  Montgomery  L&D  at  Industry, 
Pennsylvania,  is  similarto Emsworth  L&D.  The  normal 
lift  at  this  site  is  about  5.5  m.  which  promises  a  great 
potential  for  hydropower  development. 

L&D  No.  7  at  West  Kittanning.  Pennsylvania,  and 
L&D  No.  8  at  Templeton.  Pennsylvania,  on  the  Allegh¬ 
eny  River  were  visited.  Both  L&Ds  have  fixed-crest 
dams,  and  there  is  no  navigation  during  severe,  icy  win¬ 
ter  days.  Neither  heating  facilities  nor  bubblers  were 
used  there.  Only  manpower  with  chipping  hammers  and 
pike  poles  was  utilized  to  remove  solid  ice  formed  on 
the  lock  walls  and  miter  gates. 

Ice-control  measures  seem  to  depend  largely  on  each 
lockmaster's  preference.  For  example,  the  lockmaster 
at  L&D  No.  6  on  the  Mississippi  River  routinely  shakes 
all  roller  gates  to  prevent  formation  of  thick  solid  ice  on 
gate  structures.  The  frequency  of  this  loosening  opera¬ 
tion  generally  depends  on  the  weather. 

It  was  the  consensus  of  the  lockmasters  interviewed 
that,  during  typical  winters,  they  are  able  to  control  ice 
problems .  This  opinion,  again,  does  not  necessarily  rep¬ 
resent  the  general  opinion  concerning  difficulties  in 
controlling  ice  at  other  lock  and  dam  installations.  They 
all  agreed,  however,  that  recently  installed  heaters  and 
air  bubblers  had  significantly  alleviated  many  safety 
concerns.  Without  them,  manpower  was  the  only  means 
of  coping  with  ice  problems  under  extremely  treacher¬ 
ous  conditions.  Besides  safety,  operational  efficiency 
had  been  greatly  improved  with  heating  systems. 
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POWER  REQUIREMENTS 

In  this  section,  power  requirements  for  lock  wall  ice 
control  are  estimated,  and  heat-transfer  characteristics 
for  a  heat  source  embedded  in  a  concrete  lock  wall  are 
discussed. 


Heat  loss  from  a  boundary  exposed  to  frigid  air  is 
generally  estimated  in  terms  of  a  linear  combination  of 
an  ov  erall  heat-transfer  coefficient  and  the  temperature 
difference  between  the  boundary  and  air.  For  example, 
a  rate  of  heat  transfer  associated  withabandoflock  wall 
o"  j.  as  shown  in  Figure  8.  can  be  expressed  as 

°»a  =  (D 

where  //"  =  heat-loss  coefficient 

7Ws  =  wall-surface  temperature  (°C) 
fa  =  air  temperature  (°C) 

A  =  surface  area. 

A  value  of  20  \\  /in2-°C  is  commonly  used  for  the  heat- 
loss  coefficient.*  To  justify  the  use  of  this  value,  several 
references  were  surveyed.  According  to  Dickinson  and 
Chcremisinoff  ( 1 980).  heat-loss  coefficients  may  be  as 
high  as  30-60  W/m2-°C  it  an  unglazed  solar-energy 
absorber  plate  is  used  as  a  solar-energy  collector.  They 
also  reported  that  this  value  can  be  reduced  to  5-10  W/ 
m-’C  by  the  use  of  one  glazed  cover.  Therefore,  the  use 
of  20  W/m-  'C  for  the  heat-loss  coefficient  in  the 
present  case  appears  to  be  justified. 


*  Unpublished  notes.  G.D  Ashton.  1985. 


In  accordance  with  eq  I  and  with  7US  =  ()°C.  the 
powe-  . ,)ss  per  unit  area  of  lock  wall  surface,  o" ,  /A.  is 
1 00  W/m2  when  Ta  =  -5°C.  200  W/m2  when  7a  =  - 1 OC. 
and  500  W/m2  w  hen  7a  =  -25°C. 

Two  references  reporting  experiences  w  ith  heating 
of  lock  walls  to  inhibit  ice  growth  corroborate  these 
estimates  of  power  losses .  Thomas  (1976)  and  Franken¬ 
stein  and  Wortley  ( 1 984)  described  a  system  of  electric 
heating  cables  spaced  0.15  m  apart  and  embedded  in 
0.076-m  slots  cut  in  a  lock  w  all  at  St.  Lambert  Lock  on 
the  St.  Law  rence  Seaw  ay.  A  power  intensity  of  approx¬ 
imately  160  W/m2  proved  effective  in  preventing  an  ice 
collar  from  fomiing  when  the  air  temperature  was 
-5.6°C.  Balaninetal.(  1986)  reported  the  successful  use 
of  metal  heating  strips  buried  in  concrete  lock  walls. 
They  reported  that  500-600  W/m2  applied  for  three 
hours  with  an  air  temperature  of -20°C  was  usually 
sufficient  to  weaken  ice  growth  0. 1 5-0.20  m  thick  such 
that,  under  its  own  weight,  it  fell  off  lock  walls.  To 
reduce  the  power  requirements  for  de-icing.  Bulanin  et 
al.  (1986)  suggested  that  the  heating  be  performed 
intermittently  at  a  sequence  of  zones  along  the  lock 
walls. 

For  ice  control  on  lock  w  alls,  it  is  essential  tosupplv 
enough  power  to  overcome  the  power  loss  from  the  lock 
wall  boundary.  Our  study,  therefore,  set  out  to  deter¬ 
mine  the  heat  ’ransfer characteristics  of  heating  sources 
placed  in  concrete  lock  walls.  The  resulting  mathemat¬ 
ical  analysis  is  summarized  in  Appendix  A.  Several 
configurations  of  heat  sources  were  considered:  a  con¬ 
tinuous  planar  heat  source  embedded  in  the  lock  wall:  a 
line  heat  source  embedded  in  the  lock  wall:  a  rectangu¬ 
lar  boundary  containing  two  heat  sources:  and  a  triangu¬ 
lar  pocket  containing  two  line  heat  sources.  The  latter 
two  configurations  simulate  the  use  of  two  heating  rods 
inserted  in  slots  to  inhibit  ice  formation  on  roller-gate 
seals.  Each  analysis  included  temporal  distributions  of 
temperature  within  the  concrete  boundary  ,  the  temporal 
rate  of  heat  transfer  at  the  concrete-air  boundary,  and 
the  total  thermal  energy  expended  after  prescribed 
times  of  application.  The  results,  included  in  Appendix 
A.  are  presented  in  a  nomialized  format  so  as  to  be  use¬ 
ful  for  general  estimates  of  heating  requirements  for 
concrete  walls. 

A  simple,  one-dimensional,  heat-diffusion  analysis 
for  a  planar  heat  source  embedded  in  a  concrete  wall 
(Fig.  Al )  can  provide  useful  engineering  information 
on  heating  requirements  for  lock  walls.  Consider,  for 
example,  the  case  when  the  lock  wall  surface  is  at  ()°C 
and  a  heat  insulator  is  placed  at  a  =  h  (Fig.  A 1 ).  When 
a  planar  heat  source  is  placed  at  v  =  .v'.  temporal  varia¬ 
tions  of  heat  transfer  at  the  lock  wall  boundary  (v  =0)  //t 
van  be  computed  for  various  values  of  power  intensity 
(7/(6//).  in  which  Q  is  the  total  power  applied  (W).  B  is 
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Figure  9.  Temporal  variations  of  H/for 
different  Q/BH  values  (x'/b  =  0.25). 


Figure  10.  Temporal  variations  ofWffor 
different  Q/BH  values  (x'/b  =  0.50). 


T 


the  lock  wall  length  (m)  and  H  is  the  lock  wall  height 
(m).  Figure  9  shows  the  relationship  between  T  and  Hf 
for  a  range  of  QI(BH)  varying  from  50  to  300  W/m2  for 
x'/b  =  0.25,  in  which  T  is  normalized  time  [kt/(b2)],  k  is 
the  heat  diffusivity  of  concrete  (4.2  x  10 ~7  m2/s)  and  t  is 
time  (s).  Similar  relationships  are  depicted  in  Figure  10 
for  x'/b  =  0.50.  As  can  be  seen  from  these  figures,  the 
same  Hf  value  for  a  given  Q/(BH)  value  can  be  reached 
in  a  shorter  time  when  the  heat  source  is  located  closer 
to  the  lock  wall  surface. 

To  demonstrate  the  practical  use  of  these  figures,  the 
effect  of  the  power  intensity  Q/(BH)  on  the  time  re¬ 
quired  for  supplying  200  W/m2  of  heat  at  the  lock  wall 
surface  was  evaluated  for  b = 0. 1 0  m,  the  result  of  which 
is  shown  in  Table  1 .  The  relative  location  of  the  heating 
element  x'/b  was  also  varied  for  this  case.  In  accor¬ 
dance  with  eq  1 ,  Hf  =  200  W/m2  corresponds  to  a  case 
in  which  the  air  temperature  is  -10°C,  assuming  that  the 


heat-transfer  coefficient  is  20  W/m2  °C  and  Tws  =  0°C. 
Table  1  shows  that,  for  example,  it  takes  about  5.6  hours 
if  the  heater  with  an  intensity  of  300  W/m2  is  located 
C.025  m  from  the  lock  wall  surface,  and  about  7.9  hours 
if  the  intensity  is  200  W/m2.  If  the  heating  element  is 
placed  0.05  m  from  the  lock  w  all  surface  under  the  same 
conditions,  the  time  to  satisfy  H{  =  200  W/m2  is  about 
6.8  hours  and  9.2  hours  for  Ql(BH) = 300  and  200  W/m2, 
respectively.  By  using  Figures  9  and  10,  similar  esti¬ 
mates  could  be  easily  obtained  for  different  values  of  b. 

The  effect  of  the  relative  location  of  the  heating  ele¬ 
ment  on  the  rate  of  heat  transfer  at  the  lock  wall  boun¬ 
dary  was  also  analyzed.  Figure  1 1  demonstrates  the  var¬ 
iation  of  Hf  in  time  for  different  values  of  x'/b  for  Q/ 
( BH)  =  100  W/m2.  As  expected,  this  figure  shows  that 
a  larger  heat  transfer  rate  is  obtainable  with  a  smaller 
value  of  x'/b.  Table  2  illustrates  the  relative  importance 
of  x'/b  on  H(  for  b  =  0.10  m.  After  6.61  hours  of  heat 
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Figure  11.  Temporal  variations  ofH{ 
for  different  x'/b  values  (Q/BH  =  100 
W/m2). 


Table  1.  Effect  of  heat-source  intensity  and  location  of 
heating  element  on  the  time  required  for  supply  i  ng  a  rate 
of  heat  transfer  of 200  W/m2  at  the  lock  wall  surface  (one¬ 
dimensional  model). 


x'/b  =  0.25  x'/b  =  0-50 


Q/(BH) 

(W/m2) 

T 

t 

(hours) 

Q/(BH) 

(W/m2) 

T 

t 

(hours) 

50 

3.82 

25.26 

50 

4.09 

27.05 

100 

2.11 

13.96 

100 

2.35 

15.54 

200 

1.19 

7.87 

200 

1.39 

9.19 

300 

0.85 

_ 5.62... . 

..-3QQ _ 

_ LQ2 _ 

6.81 

Notes: 

Heat-transfer  rate:  Hf  =  200  W/m2. 

Plane  power  source:  Q  (W). 

Lock  wall  surface  area:  BH  (m2). 

Location  of  heating  element:  jt  =  x! . 

Concrete  lock  wall  thickness:  b  =  0.10  m. 

Heat  diffusivity  of  concrete:  k  =  4.2  x  10"7  m2/s. 
T  =  kt/(b2). 


application  at  a  power  intensity  of  1 00  W/m2,  the  rate  of 
heat  transfer  at  the  lock  wall  surface  would  be  about  94, 
8 1 , 64  and  5 1  W/m2  for  x'/b = 0. 1 0, 0.25, 0.50  and  0.75, 
respectively.  After  about  20  hours  the  rate  of  heat 
transfer  would  be  about  322, 302, 273  and  249  W/m2  for 
x'/b  =  0. 10, 0.25, 0.50  and  0.75,  respectively. 

This  analysis  is  based  on  a  one-dimensional  model 
that  utilizes  an  idealized  boundary  condition  under 
which  the  temperature  of  the  lock  wall  surface  is  as¬ 
sumed  to  be  0°C  due  to  ice  attached  to  the  lock  wall. 
Also,  it  involves  the  assumption  that  a  heat  insulator  is 
placed  at  x  =  b  such  that  the  temperature  gradient  there 
is  zero.  The  resulting  power  evaluations,  therefore,  are 
intended  for  approximately  assessing  the  required  pow¬ 
er  and  time  to  maintain  a  lock  wall  ice-free. 


Table  2.  Effect  of  heating-element  location  on  the  rate  of 
heat  transfer  at  the  lock  wall  surface  (one-dimensional 
model). 


x'/b 

Hf  (W/m2) 

t  =  6.61  hr 

1  =  13.23  hr 

t  =  19.84  hr 

t  =  26.46  hr 

0.10 

94.08 

204.18 

321.59 

443.51 

0.25 

81.37 

187.37 

302.11 

412.92 

0.50 

63.72 

162.70 

272.85 

389.11 

0.75 

51.28 

143.59 

249.08 

361.68 

Notes: 

Heat-source  strength:  Q/(BH)  =  100  W/m2. 

Lock  wall  surface  area:  BH  (m2). 

Location  of  heating  element:  x  =  x'. 

Concrete  lock  wall  thickness:  b  =  0.10  m. 

Heat  diffusivity  of  concrete:  k  =  4.2  x  IO~7  m2/s. 


FEASIBILITY  OF  UNCONVENTIONAL 
ENERGY  SOURCES 

In  the  preceding  section,  power  requirements  for  ice 
control  in  navigation  locks  were  assessed.  In  this  sec¬ 
tion,  groundwater  energy,  solar  energy,  wind  energy 
and  microhydropowerenergy  are  evaluated  as  potential 
energy  sources  for  meeting  these  power  requirements. 

Groundwater  energy 

Groundwater  energy  appears  to  be  an  attractive 
energy  source  for  use  in  ice  control  at  lock  and  dam  in¬ 
stallations.  It  is  readily  available  in  the  vicinity  of  many 
rivers,  and  year-round  its  temperature  generally  re¬ 
mains  at  the  average  annual  air  temperature  of  the  local¬ 
ity  (Fig.  12).  However,  its  attractiveness  is  somewhat 
diminished  by  practical  problems  associated  with  its 
implementation.  Additionally,  the  farther  north  (in  the 
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Figure  12.  Typical  groundwater  temperatures  in  the  United  States.  (After  Johnson  1966.) 
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Figure  13.  Typical  dimensions  of  a  lock  (Starved  Rock  Lock  and  Dam, 
Illinois  River). 


northern  hemisphere)  one  goes,  the  less  attractive  it  be¬ 
comes  because  the  groundwatertemperature  decreases. 

In  this  section,  the  power  to  heat  lock  water,  the 
power  to  heat  water  near  lock  walls,  and  the  heat- 
transfer  characteristics  of  groundwater  flow  through 
pipes  placed  in  lock  walls  are  analyzed  separately. 

Power  to  heat  water  in  locks 

Ashton*  examined  the  feasibility  ofpumping  ground- 
waterintoalockchamberwiththeobjectiveofmaintain- 
ing  the  water  temperature  above  freezing  temperature. 


*  Personal  communication.  CRREL.  1985. 


His  analysis  is  repeated  and  extended  here. 

Of  i  merest  in  the  analysis  are  the  discharge  of  ground- 
water  required  and  the  power  likely  needed  in  deliver¬ 
ing  it  to  a  representative  lock.  Ashton  assumed  that 

•  A  sufficient  quantity  of  groundwater  is  available; 

•  A  typical  groundwater  temperature  Tg  is  5.6°C, 
based  on  the  mean  annual  air  temperature; 

•  The  river  water  is  at  the  freezing  temperature,  i.e. 
Tw  =  0°C; 

•  An  appropriate  heat-transfer  coefficient  from  wa¬ 
ter  to  air  is  20  W/m2°C; 

•  Standard  lock  dimensions  are  1 83  x  33.5  m  in  plan 
and  10  m  in  depth  if  the  lock  chamber  is  held  at  the 
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Figure  14.  Heat  transfer  at  a  lock  chamber. 


upper  pool  level  or  4.5  m  in  depth  if  held  at  the 
lower  pool  level  (Fig.  13);  and 
•  No  river  water  flows  through  closed  locks. 

This  analysis  also  considers  discharge  requirements 
for  a  groundwater  temperature  of  1 2.2°C,  as  a  reason¬ 
able  maximum  value  of  groundwater  temperature  (Fig. 
1 2),  especially  for  locks  located  on  the  central  portions 
of  the  Ohio  and  Mississippi  rivers. 

The  rate  of  heat  loss  from  the  water  surface  in  the 
lock  chamber,  as  shown  in  Figure  1 4,  can  be  determined 
from 


<t>wa  =  /W7w-7a)LH/  (2) 

in  which  L  is  the  inside  length  of  the  lock  and  W  is  the 
inside  width  of  the  lock.  The  rate  of  heat  added  to  the 
lock  by  introducing  groundwater  is 

<t>g  =  Pw(?CpOrg-7'w)  (3) 

where  Q  =  discharge  of  pumped  groundwater 
pw  =  water  density 

cp  =  specific  heat  of  water  (4.2  kJ/kg  °C). 

Equating  eq  2  and  3  yields  an  expression  for  values  of 
Q  in  order  to  maintain  Tw  =  0°C  in  the  lock  under 
different  values  of  air  temperature: 

Q  =  -(H„aT3LW)/(p„cpT&).  (4) 

Table  3  lists  the  requisite  values  of  Q  associated  with  Ta 
=  -5  to  -25°C  and  Tg  =  5 .6  and  1 2.2°C.  As  expected,  less 
warm  groundwater  is  needed  for  the  larger  value  of  Tg. 

A  concern  to  be  borne  in  mind,  however,  is  that  river 
water  at  5.6°C  is  denser  than  at  0°C;  the  peak  density 
occurs  at  4°C.  Consequently  groundwater  would  sink  to 
the  bottom  of  the  lock  such  that  lock  water  would  be¬ 
come  thermally  stratified,  and  ice  may  form  over  its  sur¬ 


face.  To  prevent  this  the  lock  either  would  have  to  be 
filled  completely  with  groundwater  or  the  lock  water 
would  have  to  be  mixed  mechanically  by  means  of  air 
bubblers  or  water  jets. 

If  no  mechanical  mixing  is  provided,  groundwater 
would  cool  from  5.6°C  to  about  4°C  before  settling  to 
the  lower  portion  of  the  lock.  Under  this  condition, 
values  of  Q  in  Table  3  should  be  multiplied  by  3.5  [=5.6/ 
(5.6  -  4.0)]  to  ensure  that  Tw  at  the  water  surface  does 
not  drop  below  0°C.  If,  for  example,  Ta  =  -1 5°C,  then 
Q  =  0.275  m3/s.  At  this  discharge,  however,  a  183-m- 
long,  33.5-m-wide  and  3-m-deep  lock  could  be  filled 
with  water  at  4°C  within  a  period  of  about  18.5  hr. 
Stratification  of  water  somewhat  complicates  estimates 
of  the  requisite  discharge  to  keep  lock  water  at  or  above 
0°C.  Analytically  the  estimation  involves  a  nonlinear 
problem  insofar  as  thermal  conduction  from  a  deepen¬ 
ing,  warmer,  lower  layer  of  water  becomes  increasingly 
important.  Basic  equations  associated  with  thermal- 
energy  conservation  for  this  condition  are  presented  in 
Appendix  B.  Before  expending  effort  in  solving  this 
problem,  it  is  pertinent  to  consider  the  practicality  of 
providing  sufficient  discharge  of  groundwater.  Certain 
concerns  to  be  considered  are  as  follows: 

1 .  Magnitude  of  discharge.  As  indicated  in  Table  3, 
if  Tg  =  5.6°C,  some  0.078-0. 105  m3/s  of  groundwater 
discharge  would  be  required  when  air  temperatures 
were  -15  to  -20°C.  This  demand  is  approximately 
equivalent  to  that  for  a  community  of  10,000  people 
each  consuming  on  the  average  of  0.76  m3/day.  If  Tg  = 
1 2.2°C  the  water  demand  would  be  about  0.036-0.048 
m3/s  for  the  same  condition. 

2.  Adequacy  of  groundwater  supply.  Sufficient  aqui- 
fercapacity  is  required  to  sustain  the  requisite  discharge 
of  groundwater.  It  is  likely  that,  at  many  lock  sites, 
groundwater-recovery  wells  would  be  placed  near  a 
river  such  that  a  significant  proportion  of  groundwater 
would  indirectly  be  river  water.  A  consequence  of  this 
could  be  that  the  groundwater  temperature  may  be 
somewhat  lower  than  5.6°C  or  12.2°C.  If  groundwater 
is  not  recharged  with  river  water,  consideration  needs  to 

Table  3.  Required  groundwater  discharge  to  maintain 
lock  water  at  0°C.  (From  information  provided  by 
G.D.  Ashton  1985.) 

Air  temp  Power  loss  to  air  Groundwater  discharge  Q 
Ta  Ka  _ (mbs) _ 


ro _ <mW) _ Tg  =  5.6°C  T„  =  /2.2°C 


-5 

0.61 

0.0262 

0.0121 

-10 

1.23 

0.0523 

0.0241 

-15 

1.86 

0.0785 

0.0361 

-20 

2.45 

0.1045 

0.0480 

_=25 _ 

3.06 

_ 0.1306 _ 

_ 0.0600 _ 

Note:  1  mW  =  lO6  W  =  1340  horsepower. 
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be  given  to  the  effect  of  groundwater  removal  on  water- 
table  elevation. 

3.  Mixing  of  lock  water.  To  appreciate  the  mixing 
difficulty,  it  is  useful  to  visualize  a  discharge  of  0.078 
m3/s  into  a  lock:  picture,  say,  two  water  jets  issuing  from 
0.20-m-diameter  nozzles  at  a  mean  velocity  of  1 .3  m/s; 
each  jet  would  penetrate  about  20  m  before  the  mean  jet 
velocity  diminishes  to  about  0.15  m/s.  If  mixing  is 
assessed  in  terms  of  adensimetric  Froude  number  F&  [= 
V/(gyAp/p)l/2]  associated  with  water  flow  through  a 
lock,  with  the  mean  flow  velocity  being  (0.078  m3/s)/ 
( 10  m  x  33.5  m)  and  the  normalized  density  difference 
Ap/p  taken  as  10-4,  it  becomes  evident  that  the  lock 
water  is  not  well  mixed;  F <j  is  on  the  order  of  I0~3. 
Consequently  additional  mixing  would  be  needed  to 
ensure  that  heat  released  from  groundwater  is  rapidly 
dispersed  to  the  surface  area  of  the  lock.  If  a  manifold 
comprising  many  smaller-diameter  nozzles  (for  exam¬ 
ple,  25-mm  in  diameter)  is  considered,  concern  may 
arise  as  to  the  freezing  of  the  flowing  water  in  these 
nozzles. 

4.  Throughflow  in  the  lock.  The  groundwater  dis¬ 
charge  requirements  summarized  in  Table  3  were  ob¬ 
tained  under  the  assumption  that  river  water  does  not 
flow  through  the  lock.  However,  at  two  lock  and  dam  in¬ 
stallations  visited,  it  was  reported  that  it  is  fairly  com¬ 
mon  during  winter  to  have  lock  miter  gates  slightly  ajar, 
allowing  some  throughflow  of  water  to  inhibit  ice 
growth  within  the  central  portion  of  the  lock. 

Water  at  12.2°C  is  less  dense  than  at  0°C.  Conse¬ 
quently  it  is  more  readily  retained  at  the  surface  of  a  lock 
chamber  than  is  water  at  4°C,  at  least  until  it  is  cooled 
to  4°C.  However,  the  four  concerns  enumerated  above 
still  need  to  be  taken  into  account.  The  crucial  concern 
is  in  assuming  that  the  groundwater  becomes  uniformly 
and  rapidly  distributed  over  most  of  the  lock  chamber. 
The  power  requirements  to  pump  and  mix  groundwater 
within  a  lock  chamber  should  be  assessed  for  compar¬ 
ison  with  power  applied  for  direct  heating  of  lock  walls 
as  an  alternative  approach.  Consider  a  system  compris¬ 
ing  two  0.063-m3/s,  30.5-m  head  pumps  for  pumping 
groundwater,  and  two 0.063-m3/s,  1 5.2-m-head  pumps 
for  dispersing  groundwater  throughout  a  lock  chamber. 
The  power  requirements,  allowing  forpumpefficiency, 
would  be  about  60  kW  for  the  former  and  38  kW  for  the 
latter,  a  total  power  of  98  kW.  This  example  is  approx¬ 
imate  and  possibly  conservative,  but  it  illustrates  the 
potential  power  requirements. 

Power  to  heat  water  near  lock  walls 

An  alternative  to  heating  the  entire  body  of  water 
contained  in  a  lock  is  to  heat  water  in  the  vicinity  of  the 
lock  walls  at  the  elevation  where  ice  collars  predomi¬ 
nantly  form.  To  do  this,  a  manifold  would  be  needed  to 


Figure  15.  Heat  transfer  and  warm  water 
circulation  near  a  lock  wall. 

distribute  groundwater  flow  (or  heat)  along  the  wall, 
and  ideally  to  instigate  a  flow  circulation,  as  illustrated 
in  Figure  15.  Two  criteria  need  to  be  considered  in 
assessing  the  requisite  groundwater  discharge  to  main¬ 
tain  water  at  or  above  0°C  in  a  strip  of  width  B  along  a 
lock  wall: 

•  Thermal  power  derived  from  groundwater  dis¬ 
charge  should  balance  the  rate  of  heat  loss  to  air; 
and 

•  The  manifold  discharge  should  be  of  sufficient 
magnitude  that  manifold  jets  adequately  mix  the 
water  in  the  strip. 

The  rate  of  heat  loss  to  air  from  strips  of  width  W' 
along  both  lock  walls  can  be  expressed  as 

o;a=2/C  (Tm-TJLW.  (5) 

If  the  heat  transfer  coefficient  =  20  W/m2-°C,  W' 
=  1  m,  L  =  1 83  m  and  Tw  -  0°C,  then  eq  5  yields  values 
off)', ,  for  varying  Tv  as  summarized  in  Table  4.  These 
values  of  which  are  IW'fW  (=  2  x  1/33.5  =  0.06) 


Table  4.  Heat  loss  rate  and  groundwater  discharge 
at  strips  bordering  lock  walls. 


Air  temp 
Ta 
(°C) 

Power  loss  to  air 

Ka 

(mW) 

Groundwater  discharge  Q 
(nPls) 

T  S,=5.6V  Ti=l22‘C 

-5 

37 

0.0016 

0.0007 

-15 

no 

0.0047 

0.0021 

-25 

183 

0.0078 

0.0036 

Notes: 

1  kW  =  103  W  =  1.34  horsepower. 

Calculation  was  made  for  the  assumed  surface  area  to  be 
heated  directly:  2  x  183  m  x  I  m  =  366  m2. 
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Q/Q0  =  0.32  (h/d) 


(9) 


times  the  <}>wa  values  listed  in  Table  3,  have  to  be  bal¬ 
anced  by  heat  released  from  groundwater  discharge 
along  the  lock  walls.  The  requisite  groundwater  dis-  where  VQ  =  jet  velocity  at  the  nozzle 
charge  would  be  =  maximum  jet  velocity  at  elevation  h  above 

the  nozzle 

Q  =  <*wa  /fPwt'p(Tg  -  Tw)].  (6)  d  =  nozzle  diameter 

Bj  =  jet  width  at  elevation  h  above  the  nozzle 
Table  4  lists  values  of  Q  commensurate  with  Tg  =  5.6  Q0  =  nozzle  discharge 

and  1 2.2°C  for  a  range  of  Fa  values.  The  requisite  val-  Q  =  nozzle  discharge  plus  entrained  discharge, 

ues  of  Q  are  comparatively  modest.  They  are,  however, 

insufficient  to  mix  the  water  along  the  lock  walls.  Using  =  0.3  m/s  at,  for  example,  h  =  3  m  and  d  = 

Consider  now  the  discharge  requirements  for  a  wa-  0.025  m  yields  V0 = 5.7  m/s ,  fij  =  0.3  m,  Q0 = 2.8  x  1 0~3 

ter-jet  manifold  placed  along  a  lock  wall  so  as  to  release  m3/s  and  Q  -  0. 1 08  m3/s.  The  total  discharge  through  a 

“warm”  water  and  to  create  a  desirable  level  of  turbu-  manifold  placed  along  each  lock  wall  would  be  approxi- 

lence  at  the  water  surface.  To  estimate  the  discharge  mately  0.151  m3/s.  Evidently  much  greater  discharges 

requirements,  the  manifold  geometry  depicted  in  Fig-  are  required  for  mixing  than  for  heating.  If  mixing  is 


Figure  16.  Manifold  dimensions  for  an  assumed  diffuser  pipe. 


ure  16  is  assumed,  together  with  the  following  design 
criteria: 

•  The  maximum  velocity  at  the  water  surface,  locat¬ 
ed  at  height  h  above  manifold  nozzles,  is  0.3  m/s; 

•  The  strip  affected  by  the  circulation  cell  generated 
by  each  water  jet  is  1  m  wide;  and 

•  The  nozzles  are  0.01  m  in  diameter  and  are  spaced 
2  m  apart  along  a  1 80-m-long  manifold  pipe.  If  it  is 
assumed  that  the  total  outlet  area  should  not  exceed 
about  one  third  of  the  manifold  area,  a  0.25-  or 
0.30-m-diametermanifoldpipe  would  be  required. 

If  each  nozzle  flow  is  treated  as  a  circular  turbulent 
jet  and  the  density  difference  between  the  jet  water  and 
the  lock  water  is  neglected,  the  following  relationships 
may  be  used  to  estimate  jet  and  manifold  discharges 
(Rajaratham  1976): 


VJV0  =  6.3  (d/h) 

(7) 

fij  =  0.10/1 

(8) 

and 


assessed  in  terms  of  a  representative  densimetric  Froude 
number  F&  for  flow  through  a  lock,  it  becomes  evident 
that  lock  water  is  marginally  well  mixed;  F&  is  on  the 
order  of  unity. 

If  groundwater  at  5.6°C  is  added  at,  say,  0.0032  m3/ 
s  with  river  water  pumped  at0°C  and  0.151  m3/s  before 
being  passed  through  each  manifold,  the  blended  flow 
would  have  a  temperature  of  about  0. 1  °C.  If  the  ground- 
water  is  at  1 2.2°C,  the  blended  flow  would  have  a  tem¬ 
perature  of  about  0.25°C.  The  blended  discharge  may 
be  adequate  to  maintain  an  ice-free  strip  along  each  lock 
wall,  but  practical  difficulties  may  arise: 

•  The  blended  flow  may  cool  to  0°C  when  conveyed 
along  the  manifold; 

•  If  the  water  surface  is  located  at  an  elevation  greater 
than  3  m  above  the  manifold,  greater  water  dis¬ 
charges  would  be  needed; 

•  A  0.25-  to  0.30-m-diameter  pipe  may  itself  be  an 
unacceptable  protrusion  into  the  lock;  and 

•  The  manifold  would  have  to  be  designed  such  that 
it  can  be  quickly  drained  when  exposed  to  frigid  air 
when  the  water  level  in  the  lock  is  lowered. 
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Figure  17.  Heat  transfer  characteristics  of 
pipe  flow  (Q  =  10  gal./min). 


Heat-transfer  characteristics  of  piped  flow 
in  lock  walls 

Another  means  of  using  warm  groundwater  energy 
would  be  to  circulate  it  through  pipes  embedded  in  lock 
walls  or  through  pipes  placed  in  special  tubes  located  in 
lock  walls.  The  convective  heat  transfer  from  flow 
through  a  pipe  placed  in  a  concrete  lock  wall  can  be  esti¬ 
mated  in  the  following  manner. 

Consider  a  steady  flow  through  a  pipe  with  a  constant 
pipe- wall  temperature  T0.  The  temperature  distribution 
along  the  flow  direction  x  is  (Welty  1978) 

(r  -  r0)  /  (re  -  r0) = e~  5,(4*/D)  (io> 

where  T  =  water  temperature  at  x  (°C) 

Tt  =  entrance  water  temperature  at  x  =  0  (°C) 
St  =  Stanton  number 
x  =  longitudinal  distance  (m) 

D  =  pipe  diameter  (m). 

An  empirical  relationship  for  the  Stanton  number  devel¬ 
oped  by  Colbum  (1933)  is 

St  =  0.023  /?e~°  2  Pr~m  (11) 

where  Re  is  the  Reynolds  number  evaluated  at  the  film 
temperature,  which  is  defined  as  the  average  of  entrance 
and  exit  temperatures,  and  Pr  is  the  Prandtl  number 
evaluated  at  the  film  temperature.  7^,  the  exit  tempera¬ 
ture  at  x  =  L,  must  be  determined  from  eq  10  and  1 1  by 
a  trial-and-error  procedure  because  the  film  tempera¬ 
ture  cannot  be  predetermined.  The  exit  temperature  is 
first  assumed,  and  Re  and  Pr  are  subsequently  deter¬ 
mined.  The  Stanton  number  given  by  eq  11  is  then 
substituted  into  eq  10  to  obtain  7^.  If  this  exit  temper¬ 


ature  is  not  close  to  the  assumed  value,  the  same  pro¬ 
cedure  is  repeated  with  a  new  TL  until  the  agreement  is 
satisfactory.  Once  7"l  is  determined,  the  heat  transfer  $ 
can  be  obtained  according  to 

<|>  =  pwtAVcp(rL  -  Te))  (12) 

where  pw  =  density  of  water  (kg/m3) 

A  =  pipe  cross-section  area  (m2) 

V  =  average  flow  velocity  (m/s) 
cp  =  specific  heat  of  water  (W  s/kg  °C). 

Following  the  general  procedure  outlined  above, 
several  practical  cases  in  which  the  pipe-wall  tempera¬ 
ture  T0  was  kept  at  0qC  were  examined.  Figure  1 7  shows 
the  computed  results  of  temperature  and  heat-transfer 
variations  along  the  x  direction  in  which  D  =  0.025  m, 
Te  =  14°C,  Q-  6.3  x  10-4  m3/s  and  V=  1 .24  m/s.  As  can 
be  seen  in  Figure  1 7,  the  water  temperature  drops  very 
rapidly  along  the  flow  direction  and  reaches  the  freezing 
point  within  about  x = 200  ft.  The  heat  transferred  to  the 
lock  wall  from  warm- water  pipe  flow  was  also  estimat¬ 
ed  at  1 0-ft  intervals,  and  its  spatial  distribution  is  shown 
in  the  same  figure.  Because  the  power  <|>  is  proportional 
to  the  temperature  difference  it  decreases  with 
water  temperature.  The  total  heat  transferred  in  this  6 1  - 
m-long  pipe  segment  is  the  sum  of  <j>  shown  in  the  figure, 
which  amounts  to  about  37  kW.  Figure  18  shows  the 
computed  results  for  a  similar  case  with  a  different 
discharge  of  9.5  x  10"4  m3/s.  Figure  19  depicts  the  case 
of  Q  =  2.5  x  10-3  m3/s,  D  =  0.050  m  and  V  =  1 .24  m/s. 
The  water  temperature  decreases  rather  rapidly:  T  = 
7.3°C  at  x  =  15.2  m,  7  =  3.9°C  at*  =  30.5  m,  T  =  2.1°C 
at  *  =  45.7  m  and  T  =  1.2°C  at  *  =  61  m. 

The  total  power  expended  along  a  61 -m-long  pipe 
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T  (°C)  and  4>  (kW)  T  (°c) and  <l>  (kW)  T  (°C)  and  4>  (kW) 


Q  =  15  gpm 
— O—  X-T  Curve 
— 6---  X-<f>  Curve 

Pipe  ID.  =  I" 
T0=14°C 


Figure  18.  Heat  transfer  characteristics  of 
pipe  flow  (Q  =  15  gal. I min). 
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Q=  40  gpm 
— O—  X-T  Curve 
— 6 —  X-<f>  Curve 

Pipe  ID.  =  2“ 

T0=14«C 


Figure  19.  Heat  transfer  characteristics  of 
pipe  flow  (Q  =  40  galJmin). 
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0  =  60  gpm 
— — O—  X-T  Curve 
— Ci —  X-<f>  Curve 

Pipe  ID.  =  2" 
T0=14”C 


Figure  20.  Heat  transfer  characteristics  of 
pipe  flow  (Q  =  60  gallmin). 
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Figure  2 1 .  Distribution  of  me  an 
daily  solar  radiation  in  the 
United  States.  (After  Katiman 
1983.) 


section  is  about  136  kW.  When  the  water  discharge  is 
increased  by  50%  (2=0.0038  m?/s),  the  water  tempera¬ 
ture  at  x = 61  m  is  estimated  to  be  about  1 .4°C  (Fig.  20), 
and  the  total  power  expenditure  is  about  200  k  W.  From 
these  examples,  it  is  clear  that  unless  extremely  high 
temperature  water  (or  steam)  is  passed  through  small 
pipes  embedded  in  lock  walls,  the  water  temperature 
reaches  the  freezing  point  very  rapidly  without  travel¬ 
ing  too  far  if  the  lock  wall  temperature  is  below  the 
freezing  point.  In  conclusion,  heating  lock  walls  by 
means  of  small  pipes  carrying  groundwater  is  not  prac¬ 
tical. 

Before  determining  the  means  of  overcoming  these 
difficulties,  groundwater  use  should  be  compared  with 
other  alternative  heating  methods,  such  as  heating  ca¬ 
bles  and  steam  pipes  embedded  in  lock  walls.  Alterna¬ 
tively,  agitation  of  flow  along  lock  walls  alone  may  be 
sufficient  to  significantly  reduce  ice-collar  growth.  Riv¬ 
er  water  pumped  from  the  upper  pool  could  also  be 
passed  through  a  manifold  pipe  located  near  the  base  of 
each  lock  wall.  An  advantage  of  using  river  water  is  that 
it  is  easily  accessible  and  would  not  require  large-power 
pumps  to  circulate  it  through  a  manifold. 


Solar  energy 

Although  solar  energy  is  abundant,  even  in  the 
northern  part  of  the  U.S.,  applying  solar  energy  to  keep 
lock  and  dam  installations  ice-free  during  winter  does 
not  seem  to  be  practical .  As  can  be  seen  in  Figure  2 1 ,  ap¬ 
proximately  4  kWh/m2  of  mean  solar  energy  is  avail¬ 
able  daily  in  most  northern  midwest  and  northeast  parts 
of  the  U.S.  in  winter.  To  put  this  in  context,  consider  a 
small  house  with  a  south-facing  roof  area  of  about  1 00 
m2.  It  receives  about  400  kWh  of  solar  energy  every 
day,  or  1 2,000  kWh  per  month,  which  is  far  greater  than 
a  household’s  average  monthly  electric  use  of  1000 
kWh.  Therefore,  theoretically  speaking,  only  a  10% 
efficiency  rate  in  capturing  this  solar  energy  is  adequate 
to  meet  residential  electric  demands.  However,  this 
level  of  solar  energy  is  available  only  during  clear  day¬ 
light  conditions,  which  seldom  prevail  in  the  northern 
U.S.  because  cloudy  days  prevail  during  most  of  the 
winter.  As  a  practical  matter,  it  should  be  kept  in  mind 
that  thermal  energy  for  ice  control  is  particularly  needed 
during  nights  when  solar  radiation  is  not  available:  thus, 
air  temperature  drops  sharply. 

To  demonstrate  the  efficiency  of  solar-energy  col- 
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SOLAR  ENERGY  RECEIVED,  I^W/m2) 


Figure  22.  Efficiency  of  liquid  solar  collector. 

lectors,  a  sample  evaluation  was  made  for  a  typical, 
standard  liquid  heater  (2  x  1  x  0.15  m)  with  a  steel 
absorber,  single-glazed  glass  and  50-mm-thick  back 
insulation  (Dickinson  and  Cheremisinoff  1980).  The 
efficiency  of  a  solar  collector  e  can  be  determined  from 

e  =  FRtTa-FRHs(tl-Q/It  (13) 

where  FR  =  heat  removal  factor 

tr  =  transmissivity  (the  fraction  of  the  incom¬ 
ing  solar  radiation  that  reaches  the  absorb¬ 
ing  surface) 

a  =  absorptivity  (the  fraction  of  solar  energy 
reaching  the  surface  that  is  absorbed) 

Hs  =  overall  heat-loss  coefficient 
tj  =  liquid  inlet  temperature 
ra  =  air  temperature 

/t  =  solar  energy  received  on  the  upper  surface 
of  the  sloping  collector  surface. 

If  we  assume  that  FRt\a  =  0.73  <md  FRHS  =  4.49 


SOLAR  ENERGY  RECEIVED,  I^W/m2) 


Figure  23.  Outlet  temperature  of  liquid  solar  collector. 

W/m2-°C  (Dickinson  and  Cheremisinoff  1980),  the 
efficiency  e  can  be  estimated  using  eq  13  for  various 
combinations  of  fj,  ra  and  7t.  Table  3  summarizes  the 
computed  efficiency  for  three  air  temperatures  (-5,-10 
and  -20°C)  for  a  range  of  /t  from  300  to  500  W/m2, 
which  is  typical  around  noon  in  Minneapolis,  Minne¬ 
sota,  which  is  at  a  latitude  of  45°N.  We  assume  in  this 
estimation  that  warm  groundwater  at  10°C  is  available 
as  a  liquid  to  absorb  solar  energy.  The  results  for  effi¬ 
ciency  are  also  depicted  in  Figure  22.  The  efficiency 
decreases  very  rapidly  as  air  temperature  decreases  be¬ 
cause  heat  loss  is  directly  proportional  to  the  tempera¬ 
ture  difference  between  liquid  and  air  according  to  eq 
13.  Using  the  estimated  efficiency,  we  estimated  the 
liquid  outlet  temperature  t0  for  each  case  (Table  5,  Fig. 
23).  As  can  be  seen  in  Figure  23,  the  temperature  rise 
gained  by  means  of  aliquid  solarcollector  is  at  most  5^0 
for  a  range  of  air  temperature  from  -5  to  -20°C. 

Wind  energy 

Wind  energy  is  not  reliably  available  over  time. 


Table  5.  Computed  efficiency  and  outlet  temperature  of  liquid  solar  collector. 


Solar  radiation  l,  _ Efficiency  e _  Outlet  temperature  \0(°C) 


(W/m2) 

tfl  =  -5r 

-I0°C 

-20°C 

t„=  -5°C 

-10°C 

-20°C 

300 

0.505 

0.430 

0.280 

12.72 

12.34 

11.51 

400 

0.562 

0.506 

0.393 

13.95 

13.55 

12.82 

500 

0.595 

0.550 

0.460 

15.31 

14.92 

14.15 

Design  parameters: 

Notes: 

Fr/,  a  =  0.73. 

F* 

=  heat  removal  factor. 

FrW,  =  4.49  (W/m2- 

°C). 

t,= 

solar  radiation  transmissivity. 

Water  discharge  =  0.0 1 36  (kg/m2 

•s). 

a  = 

=  solar  radiation  absorptivity. 

Specific  heat  of  water  =  4.2  (kJ/kg-°C).  H,  =  overall  heat-loss  coefficient. 
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although  daily  and  annual  cycles  of  average  wind  speed 
are  surprisingly  stable  from  year  to  year.  The  power 
output  of  a  horizontal  windmill  strongly  depends  on 
wind  speed  (Dickinson  and  Cheremisinoff  1980): 

P  =  0.5  e  pAF3  (14) 

where  P  =  power  output  (W) 
e  =  efficiency 

p  =  air  density  =  1 .293  kg/m3  at  0°C,  1  atmo¬ 
sphere  and  sea  level 

A  =  area  of  the  rotating  components  (m2) 

V  =  wind  speed  (m/s). 

Wind  speeds  can  change  by  50%  from  one  hour  to  the 
next,  and  gusting  may  alter  wind  speeds  by  a  factor  of 
three  or  more  in  a  matter  of  seconds.  From  eq  14  it  is 
clear  that  power  varies  as  the  cube  of  the  wind  speed, 
which  means  that  a  change  in  wind  speed  by  a  factor  of 
two,  for  example,  can  result  in  a  change  in  power  output 
by  a  factor  of  eight.  It  is  not  surprising,  therefore,  that 
there  are  practical  engineering  problems  to  be  over¬ 
come  in  the  area  of  system  control,  such  that  variable 
electrical  energy  output  can  be  converted  to  constant- 
frequency  AC  that  is  phased  synchronously  with  a 
utility  network. 

Consider  the  wind-power  equation  given  by  eq  14  to 
estimate  the  power  output  that  can  be  generated  near 
Pittsburgh,  Pennsylvania,  or  Minneapolis,  Minnesota. 
In  these  areas  the  average  winter  wind  speed  is  about  4 
m/s  (Pierson  1978).  A  windmill  with  6-m-diameter 
rotor  blades  (area  =  28.27  m2)  and  an  efficiency  of  0.5 
can  generate  electric  power  output  P  as 

P  =  0.5  x  0.5  x  1.293x28.27x43  =  585  W.  (15) 

This  means  that  at  least  four  or  five  6-m-diameter 
windmill  units  are  needed  to  supply  2.5  kW  to  even  one 
unit  of  the  quartz  heaters  described  earlier.  Therefore, 
wind  energy  does  not  appear  to  be  feasible  for  supplying 
power  for  ice  control  in  navigation  locks  and  dams. 

Hydroelectric  energy 

Although  hydroelectric  energy  is  usually  thought  of 
as  a  conventional  energy  source,  it  is  included  in  this 
study  because  it  could  be  generated  in  unconventional 
ways  at  many  lock  and  dam  installations.  Microhydro¬ 
electric  power  plants  can  be  found  at  several  navigation 
lock  and  dam  installations,  such  as  L&D  No.  15  on  the 
Mississippi  River,  with  a  generator  rated  at  312  kW  at 
3.7-m  gross  head;  L&D  No.  3  on  the  Allegheny  River, 
with  a  60-kW  unit  at  4. 1-m  gross  head;  and  Emsworth 
L&D  on  the  Ohio  River,  with  a  300-kW  unit  at  5.5-m 
gross  head.  Electric  power  generated  in  these  installa¬ 


tions  has  been  used  primarily  to  operate  lock  and  dam 
gates  and  for  general  maintenance  purposes. 

Locks  and  dams  are  primarily  constructed  to  regu¬ 
late  the  stages  of  rivers  for  navigation.  The  pools  that 
they  create,  however,  can  provide  significant  hydraulic 
heads,  which  can  be  used  for  hydropower  generation. 
The  capability  to  fully  exploit  low-head  hydropower 
has,  until  recently,  been  ignored,  mainly  because  most 
available  turbines  were  developed  for  relatively  high 
heads.  However,  so-called  “low-head”  microhydropow¬ 
er  plants  have  been  built  since  the  late  1970s  at  many 
existing  low-head  dams  in  the  U.S.,  and  many  more 
projects  are  pending.  Before  discussing  the  feasibility 
of  microhydropower  plants  for  ice  control,  consider 
first  the  hydropower  output  relationship. 

The  hydroelectric  power  output  P  is  generally  ex¬ 
pressed  as 

P  =  9.%eQH  (Q  and  H  in  SI  units)  (16a) 
or 

P  =  0.0847  eQH  ( Q  and  H  in  English  units)  (16b) 

where  P  =  power  output  (kW) 
e  =  efficiency 

Q  =  water  discharge  (m3/s  in  eq  16a;  ft3/s  in  eq 
16b) 

H  =  net  head  (m  in  eq  16a;  ft  in  eq  16b). 

The  power  output,  therefore,  can  be  estimated  using  eq 
16  for  various  combinations  of  water  discharge  and 
head,  as  summarized  in  Table  6.  The  estimated  outputs 
shown  in  Table  6  were  obtained  by  assuming  that  the 
output  efficiency  is  0.8. 

Table  6  can  be  used  to  evaluate  the  feasibility  of  a 
microhydropower  plant  at  existing  lock  and  dam  instal- 

Table  6.  Hydroelectric  power  output  for  various  water 
discharge-head  combinations. 


Water  discharge,  Q 
(np/s) 

P  (kW) 

H  =  1.5  m 

3.1  m 

4.6  m 

6.1  m 

7.1 

85 

170 

255 

340 

14.2 

170 

340 

510 

680 

21.2 

255 

510 

760 

1015 

28.3 

340 

680 

1015 

1355 

34.1 

425 

845 

1270 

1695 

42.5 

510 

1015 

1525 

2035 

49.6 

595 

1185 

1780 

2370 

56.6 

680 

1355 

203 

2710 

Notes: 

P  (kW)  =  9.8  x  e  x  Q  (m3/s)  x  H  (m). 
e  =  efficiency  (assumed  to  be  0.8). 

Q  =  water  discharge. 

H  =  net  head. 
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lations.  Consider,  for  example,  L&D  No.  5  on  the 
Mississippi  River,  which  requires  about  140  kW  of 
power  to  de-ice  six  roller  gates  by  means  of  quartz 
heaters  and  side-seal  heaters,  as  described  earlier.  Be¬ 
cause  a  net  head  of  about  2. 1  m  is  available  at  L&D  No. 
5,  only  about  8.4  m3/s  of  discharge  is  sufficient  in 
generating  this  amount  of  power.  Another  example 
shows  that  if  a  discharge  of  1 4.2  m3/s  is  available  at  a  net 
head  of  3  m,  as  much  as  340  kW  of  power  output  can  be 
generated  (Table  6).  This  is  more  than  the  power  re¬ 
quired  for  heating  two  183-m-long  and  2-m-high  lock 
walls  to  control  ice-collar  growth  (approximately  300 
kW),  which  was  estimated  under  the  assumed  condi¬ 
tions  of  a  heat-loss  coefficient  of  20  W/m2oC,  a  lock- 
wall  surface  temperature  of  0°C  and  an  air  temperature 
of -20°C.  Eranti  et  al.  ( 1 983)  reported  similar  estimates 
on  the  heat  energy  required  to  maintain  satisfactory 
operation  of  lock  facilities  along  the  Saimaa  Canal  in 
Finland.  According  to  their  analysis,  approximately 
300  kW  of  thermal  energy  is  needed  for  each  lock. 

Although  unconventional  energy  sources  may  be 
cheaper  than  conventional  sources,  such  as  power  sup¬ 
plied  by  commercial  utilities,  a  crucial  question  to  be 
addressed  is:  how  can  the  requisite  level  of  hydroelec¬ 
tric  energy  be  generated  without  extensive,  and  expen¬ 
sive,  modifications  of  existing  lock  and  dam  installa¬ 
tions?  Although  the  Power  Barge™  concept  by  Broome 
(1985,  1987)  is  still  in  its  demonstration  stages,  a  pre¬ 
fabricated,  portable  power  plant  appears  to  offer  an  ex¬ 
tremely  advantageous,  and  inexpensive,  means  to  gen¬ 
erate  hydropower  for  use  in  ice  control.  There  are  two 
prototype  demonstration  projects  pending:  a  500-kW 
unit  at  a  net  head  of  5.5  m  and  discharge  of  1 1 .3  m3/s 
(Broome  1985),  anda  1 250-kW  unit  at  a  net  head  of  3.7 
m  and  discharge  of  42.5  m3/s  (Broome  1987).  The 
portable  power  plant  comprises  an  anchored  upstream 
barge  (which  provides  the  water  intake),  a  siphon  pen¬ 
stock  and  a  downstream  barge  (which  houses  a  sub¬ 
mersible  horizontal  turbine).  The  design  employs  trun¬ 
nion-type  joints  to  accommodate  variations  in  head. 
Because  no  major  construction  is  involved,  it  appears  to 
produce  minimal  environmental  impacts.  Unlike  con¬ 
ventional,  shore-attached  power  plants,  it  can  be  placed 
at  any  location  along  the  dam,  permitting  a  downstream 
flow  distribution  favorable  to  overall  navigation  and 
aquatic  environments.  This  unconventional  means  of 
generating  hydropower  energy  seems  to  be  particularly 
well  suited  for  low-head  lock  and  dam  installations. 

CONCLUSIONS 

The  principal  conclusions  derived  from  this  investi¬ 
gation  are  as  follows: 


1.  At  lock  and  dam  installations  the  primary  ice 
growths  of  concern,  and  to  be  controlled  by  heat  appli¬ 
cation,  are  ice-collar  formation  along  lock  walls,  and 
dam  gates  frozen  to  upstream  ice  sheets  and  encrusted 
by  icings  from  water  splash,  spray  and  seepage. 

2.  The  level  of  power  required  to  inhibit  ice  growth 
along  two  walls  ( 1 83  m  long  and  2  m  high)  of  a  standard 
lock  chamber  is  on  the  order  of  300  kW  when  the  air 
temperature  is  about  -20°C,  the  coefficient  of  heat 
transfer  from  lock  wall  to  frigid  air  is  taken  as  20 
W/m2-°C,  and  the  lock  walls  are  heated  to  0°C.  This 
power  requirement,  which  is  equivalent  to  about  400  W/ 
m2,  is  reduced  proportionately  with  higher  air  temper¬ 
atures.  It  can  be  reduced  further  if  lock  wall  sections  are 
heated  consecutively,  rather  than  all  of  the  walls  heated 
at  the  same  time.  Still  further  reductions  in  power  are 
possible  if  heat  is  applied  intermittently  to  dislodge  ice 
grown  on  lock  walls  rather  than  to  prevent  ice  growth. 

3.  Experience  with  typical  lock  and  dam  installa¬ 
tions,  notably  L&D  No.  5  and  6  on  the  Mississippi 
River,  indicates  that,  during  typical  winters,  roller  and 
tainter  gates  can  be  maintained  operable  by  means  of 
four  quartz  heaters  (powered  at  10  kW  total)  and  two 
side-seal  heaters  (powered  at  1 3.8  kW  total ).  Therefore, 
an  approximate  total  of  24  kW  is  required  to  keep  one 
gate  operable.  A  dam  comprising  ten  gates,  all  kept 
operable  during  winter,  would  require  about  240  k  W  for 
simultaneous  heating,  with  proportionately  less  power 
needed  if  fewer  gates  are  heated. 

4.  The  most  attractive  unconventional  energy  source 
that  meets  the  power  requirements  appears  to  be  hydro¬ 
electric  power  generated  at  locks  and  dams  by  means  of 
small,  portable,  power-barge  plants.  A  single  plant  can 
generate  300  k  W,  for  example,  at  a  net  head  of  about  3 
m  and  a  water  discharge  of  about  1 3  m3/s.  It  is  recom¬ 
mended  that  the  prototype  power-barge  projects  be 
closely  monitored  and  a  more  detailed  feasibility  study, 
including  an  economic  analysis,  be  undertaken. 

5.  Direct  heating  of  lock  walls,  lock  miter-gate  re¬ 
cesses,  and  dam-gate  pier  walls  by  means  of  electrical 
heaters  embedded  in  concrete  walls  is  the  best  approach 
for  controlling  ice  growth  on  lock  walls  and  dam  gates. 

6.  Neither  solar  energy  nor  wind  energy  is  feasible 
for  ice  control  at  lock  and  dam  installations  because 
they  are  unreliable,  since  no  solar  power  is  generated 
during  nights  and  no  wind  power  is  generated  during 
windless  days. 

7.  The  use  of  groundwater  tocontrol  ice  formation  in 
locks  seems  attractive.  However,  the  practical  aspects 
of  mixing  groundwater  with  lock  water,  or  discharging 
it  along  lock  walls,  diminish  its  general  utility  for  ice 
control.  Additionally,  as  the  groundwater  temperature 
decreases  with  increasing  latitude,  so  does  the  attrac¬ 
tiveness  of  its  use  for  ice  control.  The  analysis  of  heating 
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water  near  lock  walls  also  shows  that  a  greater  amount 
of  discharge  would  be  required  for  mixing  than  for  heat¬ 
ing  water  along  each  lock  wall.  It  is  not  feasible  to  use 
warm  groundwater  by  passing  it  through  pipes  em¬ 
bedded  in  lock  walls  for  maintaining  lock  walls  ice-free. 
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APPENDIX  A:  MATHEMATICAL  ANALYSIS  OF 
HEAT-TRANSFER  PROCESSES  THROUGH  NAVIGATION  LOCK  WALLS 


ONE-DIMENSIONAL  HEAT-DIFFUSION  MODEL 


Because  the  heat  conductivity  of  concrete  is  extremely  low  (approximately  two  orders  of 
magnitude  smaller  than  that  of  metal),  the  use  of  an  isolated  heating  rod  in  a  two-dimensional 
field  as  a  point  source  is  not  effective  in  heating  a  long  lock  wall.  Therefore,  a  one¬ 
dimensional  heat-transfer  model  with  a  continuous  plane  heat  source  along  the  lock  wall  was 
first  evaluated. 

Consider  a  continuous,  steady  heat  source  embedded  in  concrete  along  a  lock  wall,  as 
shown  in  Figure  A1 .  The  governing  equation  for  a  one-dimensional  heat  diffusion  is 

?!L  =  k—  (Al) 

to  dx2 

where  0  =  temperature  (°C) 

k  =  heat  diffusivity  (m2/s) 
x  =  distance  (m) 
i  =  time  (s). 


If  heat  is  liberated  at  the  rate  of  p cpq  (W/m2)  per  unit  area  per  unit  time  in  the  plane  x  =  x' 
at  time  t  =  0,  a  general  solution  of  eq  Al  is 


0  =  q  |(f/jtk) 


1/2 


exp 


(*-*f 

4kt  J 


^'Urfc^-*' 


2k 


2  VF 


=  E{x- x') 


(A2) 


Figure  Al .  Definition  sketch  of  plane  heat-source  model. 
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where  q  =  heat-source  strength  (°Cm/s) 
p  =  density  of  solid  (kg/m3) 

C  p  =  specific  heat  of  solid  (W  s/kg-°C) 

E  =  function  expressing  temperature  distribution  when  the  heat  source  is  placed  at  .v 
/ 

=  A  . 


If  the  temperature  of  the  lock  wall  surface  is  assumed  to  be  0°C  due  to  ice  attached  to  the  lock 
wall,  and  a  heat  insulator  is  placed  at  .v  =  b,  as  shown  in  Figure  Al,  the  boundary  conditions 
become 


0  =  O°C 


at  x  =  0  (/  >  0) 


and 


^  =  0  atx-  b(t>  0). 

3a 


(A3) 


To  satisfy  eq  3  it  is  necessary  to  superimpose  a  plane  sink  at  x  =  -x'  and  a  plane  source  at  (,v 
=  2)  (b-x'),  which  yields  the  following  temperature  distribution: 


0  =  £(.v  -  a')  -  E( x  +  a')  +  E( a  +  a  -  2b).  (A4) 

It  is  more  convenient,  for  the  purpose  of  generalizing  the  analysis,  to  use  the  following 
normalized  variables: 


X  =  x/b 


X'  =  x'tb  • 
T  =  ktlb2. 


Normalized  temperature  0,  corresponding  to  eq  A2,  can  be  written  as 


0  =  —5—  =  f  T-  exp 

q(b/k)  v  n 
Similarly,  eq  A4  becomes 


0  =  — ^ — ■  =  v,  T-  exp 


AT  J  2  \  iff 


q(b/k) 


(X-X'f 


AT 


Lv-rl-erfc(i^l 
2  ff 


(A5) 


(A2a) 


+ 


(x+r)2 

+  ^X+X' 1  erfc 

||X+X'|| 

AT 

2 

iff  / 

(X+X'-2f 

|X+X'-2|erfc 

|X+X'-2  I  \ 

AT 

2 

iff  / 

(A4a) 


Temporal  variations  of  0,  as  predicted  using  eq  A4a,  are  shown  in  Figures  A2  and  A3  for  two 
plane-source  positions:  X'  =  0.25  and  0.50,  respectively.  The  information  in  Figures  A2  and 
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Figure  A2.  Spatial  distributions  of  normalized  temperature  in  one -dimen¬ 
sional  model  (x'/b  =  0.25). 


X 


Figure  A3.  Spatial  distributions  of  normalized  temperature  in  one-dimen¬ 
sional  model  (x'/b  =  0.5) 


A3  can  be  used  to  estimate  power  requirements  for  keeping  a  lock  wall  ice-free.  Taking  the 
height  of  the  zone  of  ice  formation  on  a  lock  wall  of  length  B  to  be  H,  and  the  power  supplied 
to  the  plane  source  Q,  the  requisite  heat  strength  q  to  maintain  ice-free  conditions  over  the 
wall  is 


and 


q= 


Q 

PCpBH 


(A5) 


Qb  \ 

KBH< 

where  K  is  the  thermal  conductivity  (pcp£)  in  W/m  °C. 

Consider  a  lock  wall  under  the  following  conditions:  B  =  250  m,H  =  2m,b  =  0.30  m,  k 
(concrete  wall) = 4.2  x  1 0-7  m2/s,  a  '=0.075  m  and  Q = 250 kW.  The  approximate  temperature 
at  x  =  0.03  m  when  t  =  2.4  hours  can  be  obtained  as  follows: 


Normalized  heat-source  location:  X'  =  x lb  =  0.075/0.30  =  0.25 
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Normalized  time: 


Normalized  location: 


j-kt  _  4.2  x  IQ-7  x  2.4  x  3600  _  q  ^ 
b 2  0.3  2 

X=x/b  =  0.03/0.30  =  0.10 


Normalized  temperature:  0  =  0.038  (from  Figure  A2). 

Thus,  the  actual  temperature  0  is 

e_6feg_  0.038  [0.03  (m)]  [250  x  103  (W)]  _670r. 
KBH  [0.85  (W/m  •  °C>]  [250  (m)]  [2  (m)] 


The  power  per  unit  area  of  wall  //f  associated  with  heat  transferred  from  the  concrete 
surface  to  balance  ice  formation  at  the  lock  wall  (x  =  0)  is  then 

H(=  -K  —  (w/m2) .  (A7) 

dx 

Therefore,  the  total  heat  transfer  Ha  in  time  t  can  be  determined  as 


//ft  =  |  |-AT  dr  (w-hr/m2) 


(A8) 


which  can  be  normalized  as 


(A9) 


Figure  A4  shows  Hfl  plotted  as  a  function  of  T  for  the  case  of  a  heating  plate  placed  in  the 


-i 


Figure  A4.  Temporal  variation  ofH ft  in  one-dimensional  model  (x'/b  =  05). 
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middle  of  the  lock  wall  ( X '  =  0.5).  The 
cumulative  heat-transfer  quantity  //ft  un¬ 
der  this  condition  increases  monotonical- 
ly  with  T. 

Table  A1  lists  temporal  variations  of 
the  total  heat  flux  //ft,  which  were  com¬ 
puted  for  several  combinations  of  x'  and 
BH  (lock  wall  surface  area  protected  from 
ice  formation)  when  Q  =  250  kW.  To 
estimate  the  rate  of  heat  transfer  required 
to  keep  the  concrete  wall  surface  ice-free, 
the  temperature  gradients  that  appear  in 
eq  A7  were  evaluated  at  x  =  0.  Figure  AS 
shows  temporal  variations  of  the  normal¬ 
ized  temperature  gradient  for  three  values 
of  x'/b.  As  can  be  seen  in  this  figure,  the 


Table  Al.  Summary  of  computed  total  amount  of 
heat  transferred  at  *  =  0,  ffn  (W-hr/m2). 

Time 


(hours) 

BH  =  500  mi2 

BH  =  250  n?  BH  =  SO  n? 

Case  UQ 

s 

© 

a 

ii 

b-- 

■  0.3  m;  X'  = 

0.15  m 

1 

0.6 

1.2 

5.8 

2 

11.6 

23.3 

116.3 

3 

58.1 

116.3 

581.4 

4 

116.3 

232.6 

1162.8 

5 

232.6 

465.1 

2325.6 

Case  2:  Q 

=  250  kW; 

b  = 

;  0.2  m;  x’  = 

0.10  m 

1 

9.3 

18.6 

93.0 

2 

79.1 

204.7 

790.7 

3 

158.1 

317.4 

1581.4 

4 

303.5 

608.1 

3034.9 

5 

594.2 

1188.4 

5941.9 

closer  the  heat  source  is  to  the  wall  surface,  the  larger  the  temperature  gradient  is  at  a  given 


time.  The  heat  flux  Hf  can  then  be  calculated  by 


H{  =  -K —  =  . 

dx  BH  dX 


(A10) 


Relationships  between  T  and 
Hf  computed  for  five  values 
for  QIBH  using  eq  A10  are 
shown  in  Figure  A6.  These 
curves  were  obtained  for  the 
case  of  a  heating  plate  placed 
in  the  middle  of  the  wall  (x'/b 
=  0.5). 


Figure  A5.  Temporal  varia¬ 
tions  of  temperature  gradient 
for  different  x'/b  values. 


Figure  A6.  Temporal  varia¬ 
tion  ofHffor  different  Q/BH 
values  (x'/b  =  0.5). 
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The  simple  one-dimensional, 
plane-source  model  presented  here 
was  investigated  to  demonstrate 
its  capability  for  estimating  the 
necessary  amount  of  heat  and  time 
of  application  that  are  required  to 
maintain  the  lock  wall  ice-free. 

Needless  to  say,  this  model  is  crude 
and  makes  various  assumptions  to 
simplify  the  problem  to  obtain 
practical,  simple  engineering  so¬ 
lutions. 

TWO-DIMENSIONAL 
HEAT-DIFFUSION  MODEL 

Consider  a  two-dimensional 
temperature  distribution  for  which 
a  steady  line  heat  source  is  placed 
in  a  concrete  lock  wall  parallel  to 

the x-axis  direction  (perpendicular  to  the  x-y  plane),  as  shown  in  Figure  A7.  It  is  assumed  that 
the  line  source  is  located  at  a  point  (a,0)  and  that  the  continuous  steady  heat  is  liberated  at  a 
rate  of  pcpq  per  unit  length  and  per  unit  time  (W/m).  The  governing  equation  in  this  case  is 


Figure  A7.  Definition  sketch  of  line-heat-source  model. 


ae  + 

3f  la*2  dy2l 


(All) 


A  6„.r-ral  solution  of  eq  A1 1  is 


e  =  -JL 

4nk 


l  e^du=EiU 1) 

J  jrM  u  Uktl 

\4kt  / 


(A12) 


where  6  =  temperature  (°C) 

q  =  line-source  strength  (°C  m2/s) 
k  -  heat  diffusivity  (m2/s) 
t  =  time  (s) 
u  =  dummy  variable 
r2  =  (x-a)2  +  y2 
E\  =  exponential  integral. 


If  the  same  boundary  conditions  (eq  A3)  as  used  in  the  previous  case  are  imposed,  it  is 
necessary  to  superimpose  a  line  sink  at  (-a,0)  to  satisfy  the  boundary  condition  at  x = 0,  and 
a  line  sink  and  a  source  at  ( 2b+a,0 )  and  (2b-a,0),  respectively,  to  satisfy  the  boundary 
condition  at  x  =  b.  The  solution  to  eq  A1 2  is 


0  =  E: 


(x-af+y: 


4kt 


-  E: 


(x+af+y: 


4  kt 


+  E; 


(. x-2b+af+y 2 


4  kt 


-  E: 


{x-2b-af+y2 


4kt 


(A13) 


26 


This  equation  can  be  normalized  as 


e=(ip 

Unkl 


{X-a/bf+Y: 


4  T 


-E{ 


[X+a/bf+Y2]  \{X-2 +a/bf+Y 


4  T 


4T 


{X-2-albf+Y2 
4  T 


(A  14) 


where  X  =  xlb 
Y  =  ylb 
T  =  ktlb2. 

Distributions  of  the  normalized  temperature  0  at  T=  0. 1  are  shown  in  Figure  A8  for  alb = 0.5 
(the  line  source  is  located  in  the  middle  of  the  lock  wall)  and  three  values  of  Y:  0, 0.25  and 
0.50.  Because  of  the  low  heat  conductivity  of  concrete,  heat  transfer  is  extremely  poor  in  the 
vertical  (y)  direction.  For  example,  at  Y  =  0.5  (horizontal  plane  that  is  0.5 b  above  the  axis), 
9  is  less  than  one-third  of  its  value  at  X  =  0. 1 5  in  comparison  with  that  when  K = 0  (Fig.  A7). 
The  value  of  9  for  Y  =  0  at  X  =  0.5  is  approximately  3 1 .78,  and  the  curve  in  the  vicinity  of  X 
=  0.5  is  not  shown  in  Figure  A8.  Figure  A9  shows  lateral  distributions  of  0  at  T  =  0.10  and 


X 


Figure  A8.  Spatial  distributions  of  normalized  temperature  in  two- 
dimensional  model  (a/b  =  0.5  and  T  =  0.10). 


Figure  A9.  Spatial  distributions  of  normalized  temperature  in  two-dimen¬ 
sional  model  (T  =  0.10  and  Y  =  0). 
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Figure  A10.  Temporal  distributions  of  //ft  (Y  =  0). 


Figure  A 12.  Temporal  distributions  of  (a/b  =  0.50). 


Figure  Al  l  .Temporal  distributions  ofHh  (a/b  =  0.25). 

K=  0  for  three  values  of  alb :  =  0.25, 0.50  and  0.75.  The 
line  heat  source  located  closer  to  the  lock  wall  (see  the 
curve  for  alb  =  0.25,  for  example)  produces  a  larger 
temperature  gradient  at  X = 0  than  those  located  farther 
from  the  lock  wall.  This  means  that  more  heat  flux  is 
available  to  melt  ice  at  a  given  time  with  a  smaller  value 
of  alb.  At  T  =  0. 1 ,  values  of  0  near  the  face  of  the  lock 
wall  (X=  0. 1 )  are  1 .467, 0.444  and  0. 145  for  alb  =  0.25, 
0.50  and  0.75  respectively.  Heat  tends  to  build  up  at  X 
=  1  for  a  larger  value  of  alb. 

To  estimate  the  rate  of  heat  transfer  from  the  concrete 
lock  wall  to  ice,  the  heat  flux  Hf  atX  =  0  was  calculated 
using  the  following  expression: 

H{=  -K  —  (A15) 

dx  x=0 


where  H(  is  in  W/m2  and  K  is  the  thermal  conductivity 
in  W/m-°C.  The  phase  change  of  ice  due  to  conducted 
heat  at  the  lock  wall  boundary  was  not  considered  for 
simplicity’s  sake.  The  total  heat  flux  //ft  accumulated  in 
time  /  at  X  =  0  can  be  obtained  by  integrating  H{  as 
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where  H ft  is  in  W  hr/m2.  Introducing  eq  A13  into  eq  A16  yields,  after  some  mathematical 
manipulations,  the  following  solution  for  //ft : 

Pcp<kb 


2n 


exp 


A2+Y2 
4  T 


,  (M 

(2-Af+Y2 


exp 


(2^) 

(2+Af+Y2 


exp 


(2-A  f+Y  2 
4  T 


(2+Af+Y2 
4  T 


df 


(A17) 


where  A  =  alb.  Figure  A 1 0  presents  the  relationships  between  //ft  and  T  at  Y = 0  for  alb = 0.25 
and  0.50,  and  shows  the  sensitivity  of  Hft  to  the  location  of  a  heat  source.  By  placing  a  line 
source  at  the  quarter  point  from  the  lock  wall  (alb = 0.25),  the  total  heat  flux  Hft  accumulated 
in  time  T  is  seen  to  be  approximately  ten  times  greater  than  that  for  the  case  of  alb  =  0.50. 
Figure  A1 1  shows  the  temporal  variations  of  H{x  for  alb  =  0.25  at  three  vertical  locations  (Y 
=  0, 0.5  and  1 .()).  The  total  heat  flux  is  seen  to  decrease  drastically  in  the  vertical  direction; 
for  example,  H{1  at  F  =  y/b  =  0.5  (along  a  line  on  the  lock  wall,  half  of  the  concrete  wall 
thickness  b  above  or  below  the  line  source)  is  almost  ten  times  smaller  than  that  at  Y = 0  and 
more  than  hundred  times  smaller  at  Y  =  1 .0  for  a  range  of  7"  between  0  and  0.20.  Similar  trends 
are  seen  in  Figure  A12  when  the  line  source  is  located  along  the  center  of  the  lock  wall  (al 
b  =  0.5).  This  analysis  clearly  demonstrates  that  the  two-dimensional  line  heat  source  is  not 
effective  in  maintaining  ice-free  conditions  over  a  significant  portion  of  the  lock  wall  because 
the  majority  of  the  heat  transfer  at  the  lock  wall  boundary  takes  place  only  in  the  vicinity  of 
Y  =  0.  Table  A2  summarizes  the  results  of  sample  computations  of  the  total  heat  flux  Hft  (W 
hr/m2)  with  a  line-source  strength  of  q  =  1 1 ,8°C  m2/s,  which  corresponds  to,  for  example,  a 
case  in  which  Q  =  250  kW  is  applied  to  a  1 00-m-long  line  source.  As  stated  previously,  there 
can  be  seen  strong  gradients  of  the  heat  flux  in  the  vertical  direction  in  this  table. 


Table  A2.  Summary  of  computed  total  amount  of  heat  transfer  Hr 
(W  hr/m2)  for  a  line  heat-source  model. 


Time 

CASE  l  :  a  = 

:  0.15  m.b  =  0.30  m 

CASE  2:  a  = 

0.075  m.b  =  0.30  m 

(hours) 

y  =  0m 

v  =  0.15  m 

v  =0  m 

y  =  0.15  m 

2 

442 

23 

6,395 

105 

4 

3.721 

500 

23.256 

1,047 

6 

9,419 

1,628 

45,349 

2,674 

Time 

CASE  3:  a  = 

>  0.10  m.b  -  020  m 

CASE  4:  a  = 

=  0.05  m,  b  =  0 20  m 

(hours) 

y  =  0  m 

y  =  0.10  m 

y  =  0  m 

y  =  0.10  m 

1 

279 

23 

4,419 

81 

2 

2,209 

291 

12,791 

640 

3 

4,535 

965 

20,698 

1.512 

4 

8.372 

2.093 

28,605 

2.674 

5 

11,628 

3,140 

41,860 

_ MM _ 

Note:  q  =  1  l.8°C  m2/s. 
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TWO-DIMENSIONAL 
HEAT-DIFFUSION  MODEL  IN 
A  RECTANGULAR  BOUNDARY 


To  heat  roller-gate  seals,  as  used  for  most 
lock  and  dam  installations  along  the  upper 
Mississippi  River,  it  is  common  to  provide 
concrete  support  piers  with  a  vertical  rectan¬ 
gular  pocket  in  which  two  heating  rods  are 
inserted  to  inhibit  ice  formation  on  the  gate 
seals.  A  typical  configuration  is  shown  in  Fig¬ 
ure  A 1 3 .  The  rectangular  pocket  is  assumed  to 
be  surrounded  by  a  metal  conductor  at  the  lock 
wall  surface  and  three  insulated  walls. 

The  governing  equation  for  heat  flux  from 
the  heating  rods  is  given  in  eq  A 1 1 ;  however, 
new  boundary  conditions  must  be  imposed  as 
follows: 


10,  d) 


•* 

— INSULATOR 

_ heat  SOURCE 

B 

m 

(o,c) 

X 

pH 

i 

.0/ 

-  •  -  », 

- CONCRETE 

-  HEAT  SOURCE 

CONDUCTOR 

m 

•  •  . 

.  ••  • 

m 

10, -d) 

m 

b 

Figure  A13.  Definition  sketch  of  line-heat- 
source  model  in  a  rectangular  boundary. 


9  =  0°C  at  a  =  0;  <  y  <  d  (/  >  0) 


^.  =  0 
9a 


at  x  =b;  -d<  y  <  d(t  <  0) 


and 


at  0  5  a  <  b\  y-±d(t>  0). 


(A  18) 


If  two  line  heat  sources  are  placed  at  (a,c)  and  (a,~c),  Sommerfield’s  method  of  images  (Carslaw 
and  Jaeger  1947)  can  be  applied,  as  shown  in  Figure  A14,  to  determine  the  temperature 
distribution  throughout  the  heated  pocket.  An  expression  for  normalized  heat  flux  H{  at  X  =  0 
can  be  obtained  in  a  similar  manner  to  that  used  earlier.  The  expression  is 
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W+Y\ 


A  2+Y?  | 
4  T  ' 


(2-A) 

(2  -Af+Y? 

-  (2+/1)  e-D 

(2+A)2+y2l| 

(2-A)2+Fi2 

i  4  T  j 

(2+A)2+yj2 

AT  Jj 

(A  19) 


where  Y 
T 

Q 

i 

b 

H 

Hi 

A 

Yi 

Yl 

Yl 

Y4 


y/b 
kt/b 2 

total  heat-source  strength  (W) 

subscript  for  location  of  image 

depth  of  the  rectangular  pocket  (m) 

height  of  ice-covered  portion  of  lock  wall  (m) 

heat  flux  at  a  =  0  (W/m2) 

alb 

Y-clb 

Y-2d/b  +  c/b 

Y  +  clb 

Y  +  2d/b  -  c/b. 
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y 


(-a,-c  +  2d) 


(o,-c+2d)  (-o  +  2b,-c+2d)  (o+2b,-c  +  2d) 


r- 

+ 

n 

' 

. 

+ 

4- 

(-o,c) 

(o,c) 

(-a  +  2b,c) 

(o  +  2b,  c) 

(-o,-c) 

(a,-c) 

<-a  +  2b,-c) 

(a  +  2  b,  -c) 

' 

+ 

+ 

_ 

+ 

+ 

(-o,c- 

2d) 

(a,c-2d)  (-a+2b,c-2d) 

(o  +  2b,  c-  2d) 

Figure  A14.  Locations  of  imaginary  sinks 
and  sources for  the  line-heat-source  mod¬ 
el  in  a  rectangular  boundary. 


The  normalized  total  amount  of  heat  expended  Hh  is  determined  by  integrating  eq  A 1 9  over  T: 


"f iMs 


dr. 


(A20) 


The  area-averaged  normalized  total  amount  of  heat  expended  Hfla  is 


(A21) 


//f(a  was  numerically  integrated,  the  results  of  which  are  presented  in  Figures  A15  and  A16  as 
a  function  of  T.  Figure  A 1 5  shows  the  effect  of  dlb  (the  aspect  ratio  of  the  rectangular  pocket) 
for  fixed  values  of  alb  =  0.5  (heat  sources  located  along  the  centerline  of  the  pocket)  and  c/b  = 
2.0.  The  aspect  ratio  is  important  for  very  small  values  of  T;  however,  the  effect  is  seen  in  the 
figure  to  diminish  for  T  greater  than  about  0.05.  Figure  A 1 6  similarly  shows  the  effect  of  c/b  (the 
positions  of  the  rods  in  the  y  direction)  for  alb  =  0.5  and  dlb  =  4.0.  These  results  demonstrate 
that  //fta  is  not  significantly  influenced,  except  for  the  time  period  close  to  heat  application,  by 
either  the  location  of  the  rods  or  the  length  of  the  pocket  for  the  given  pocket  depth  b. 

Table  A3  depicts  a  practical  example  of  how  much  heat  is  transferred  through  the  metal 
conductor  of  the  pocket  in  the  lock  wall  under  the  following  conditions:  Q-  5  kW,  H  =  2  m, 
a  =  0.05  m,  b  =  0. 1 0  m,  c  =  0.20  m  and  d  =  0.40  m.  The  total  amount  of  heat  expended  is 
shown  in  the  table  as  a  function  of  time  t.  According  to  this  table,  approximately  1.18  kW  hr/ 
m2  of  heat  can  be  accumulated  in  approximately  60  minutes.  According  to  the  Iockmaster  at 
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0,20  Figure  A16.  Temporal  distributions 
of  Hfa  fa/b  =  0.50  and  d/b  =  4.0). 


Lock  and  Dam  6  located  in  Trempealeau, 
Wisconsin,  along  the  Mississippi  River,  the 
de-icing  task  with  side-seal  heaters  can  be 
accomplished  within  about  45  minutes  of 
their  application.  Therefore,  the  estimation 
presented  above  is  quite  encouraging,  al¬ 
though  the  analytical  solution  was  obtained 
under  various  simplifying  assumptions.  Also 
encouraging  is  a  study  that  supports  the 
present  estimate,  that  of  Balanin  et  al.  ( 1 986), 
who  found  that  it  took  about  3  hours  under  - 
20°C  air  temperature  to  remove  thick  ice 
(0. 1 5-0.20  m  thick)  from  the  wall  by  apply¬ 
ing  1 . 5-2.0  kW  hr/m2  of  electricity  by  means 
of  metal  strips  buried  in  the  concrete  wall. 
According  to  Balanin  et  al.  (1986),  thick  ice 
formed  on  the  lock  wall  tends  to  slide  off  the 
wall  surface  due  to  gravity  even  before  the 
temperature  of  the  concrete  wall  surface 
reaches  0°C.  The  results  of  the  present  study 
shown  in  Table  A3  are  practically  on  the 
same  order  of  magnitude  in  all  variables  in¬ 
volved.  Therefore,  the  simple  analytical  so¬ 
lution  developed  in  the  present  study  could 
be  used  effectively  in  estimating  the  heating 
efficiency  of  the  method  of  heating  side  seals 
commonly  practiced  by  the  U.S.  Army  Corps 
of  Engineers. 

TWO-DIMENSIONAL 
HEAT-DIFFUSION  MODEL 
IN  A  TRIANGULAR  BOUNDARY 


Table  A3.  Summary  of  computed  total 
amount  of  heat  transfer  for  a  two-dimen¬ 
sional  model  in  rectangular  boundary. 


Time,  t  Dimensionless 

H 'f,o 

(min) 

lime.  T 

H f,a 

(W-hr/m2) 

10 

0.025 

0.00032 

6.6 

20 

0.050 

0.0040 

82.6 

30 

0.075 

0.012 

247.7 

40 

0.100 

0.023 

474.9 

50 

0.125 

0.039 

805.2 

60 

0.150 

0.057 

1,176.7 

a  =  0.05  m 

H  = 

2  m 

b  =  0.10  m 

*  = 

4.2  x  10-7 

m2/s 

c  =  0.20  m 

e  = 

5  kW 

d  =  0.40  m 

Figure  A17.  Definition  sketch  of  line-heat  - 
source  model  in  a  triangular  boundary. 


Consider  two  line  heat  sources  that  are  placed  in  a  triangular  domain  of  a  lock  wall  in  which 
the  hypotenuse  of  an  isosceles  right  triangle  is  the  ice-concrete  boundary  made  of  a  metal 
conductor,  as  shown  in  Figure  A 17.  The  governing  equation  is  given  by  eq  Al  1,  and  the 
boundary  conditions  are 


0  =  °C  ata:  =  0;-62y<b(r>0) 
and 


(A22) 


—  =  0  along  insulated  boundary  J 
dn 

where  n  is  normal  to  the  sides  of  the  triangle.  To  satisfy  these  boundary  conditions,  imaginary 
sources  and  sinks  are  necessary,  as  shown  in  Figure  A 18.  The  normalized  heat  flux  Ht  is 


HlY,T)  =  -?- 1 


=  -±{1 


IQ]  2n  \i=l  e*+{Y-mf 
\bH' 


exp 


4T 


(A23) 
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where  Y  =  y/b 
T  =  kt/b 2 

i  =  subscript  for  location  of  the  image 
h  =  half  of  the  wall  width 
It  =  height  of  the  wall 
£i  =  i2  =  alb 
f  3  =  (4  =  l  -  db 
—  ((,  =  \  +  db 
mi  =  -  m2  -  db 
W3  =  -  W4  =  1  -  alb 
m$  -  -ni(,  =  alb  -  1 . 


Therefore,  the  normalized  total  amount  of  heat  H(t 
expended  over  normalized  time  T  is 


T 

H^Y,T)dT 


0 


(A24) 


The  area-averaged  normalized  total  amount  of 
heat  Wfta  is 


(A25) 


H fta  was  computed  as  a  function  of  T  for  three 
cases,  and  the  results  are  shown  in  Figure  A 19. 
Case  3  in  Figure  A 1 9  represents  the  case  in  which 
a  single  rod  is  placed  at  (a,0)  where  a  -  0.5 b.  The 
single  heat  source  is  less  efficient  than  those  cases 
with  two  heating  rods,  in  particular,  for  small  T. 
Table  A4  lists  the  results  of  the  computed  total  heat 


Table  A4.  Summary  of  computed  total 
amount  of  heat  transfer  for  a  two- 
dimensional  model  in  triangular  bound¬ 
ary. 


T 

l 

(mini 

Hfta 

tyira 

(W  hr/m2) 

Case  1:  alb 

=  0.25,  db 

=  0.25 

0.010 

15.9 

9.8E-5 

16.2 

0.015 

27.8 

4.3E-4 

71.0 

0.020 

31.7 

1.0E-3 

165.2 

0.025 

39.7 

1.9E-3 

313.8 

0.030 

47.6 

3.0E-3 

495.6 

0.035 

55.6 

4.2D-3 

693.8 

0.040 

63.5 

5.7E-3 

941.6 

0.050 

79.4 

9.0E-3 

1486.7 

Case  2:  alb 

=  0  JO.  db 

=  0.30 

0.010 

15.9 

I.8E-5 

3.0 

0.015 

23.8 

I.2E-4 

19.8 

0.020 

31.7 

3.6E-4 

59.4 

0.025 

39.7 

7.6E-4 

125.6 

0.030 

47.6 

1.3E-3 

214.8 

0.035 

55.6 

2.0E-3 

330.3 

0.040 

63  s 

2.9E-3 

479.1 

0.050 

79.4 

5.0E-3 

825.9 

Case  3:  alb 

=  0  JO,  db 

=  0 

0.010 

15.9 

3.5E-6 

0.6 

0.015 

23.8 

3.9E-5 

6.4 

0.020 

31.7 

I.5E-4 

24.8 

0.025 

39.7 

3.8E-4 

62.8 

0.030 

47.6 

7.5E-4 

123.8 

0.035 

55.6 

1.2E-3 

198.3 

0.040 

63.5 

1. 9E-3 

313.8 

0.050 

79.4 

3.5E-3 

578.1 

b  =  0.20  m 
//  =  2  m 
Q=  5kW 
*  =  4.2  x  I  O'7  m’/s. 


(-b-c,  b-o) 

1  y 

(0,  b) 

n.  (b-c,  b-o)  (b  +  c,b-o) 

(-a,c) 

BIB 

1 

0 

(-a 

.-c) 

(-b-c,  a-b) 

y^b-c.o-b)  (b+c.a-b) 

(0,-b) 

Figure  A18.  Locations 
of  imaginary  sinks  and 
sources  for  the  line- 
heat-source  model  in  a 
triangular  boundary. 
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Figure  A19.  Temporal  distribu¬ 
tions  ofH^Jor  different  a/b  and 
c/b  values. 


T 


transferred  in  time  for  each  of  the  three  cases  presented  in  Figure  A 19.  It  was  assumed  that 
b  =  0.20  m,  h  =  2  m  and  (2  =  5  kW  (or  2.5  kW  per  each  heating  rod).  For  example,  the  total 
amounts  of  heat  transferred  at  the  wall  boundary  after  applying  heat  for  40  minutes  ( T  = 
0.025)  are  about  314,  126  and  63  Whr/m^  for  cases  l,2and  3,  respectively. 
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APPENDIX  B:  MATHEMATICAL  FORMULATION  OF 
INCOMPLETE  MIXING  OF  GROUNDWATER  IN  A  LOCK  CHAMBER 


The  system  of  heat  fluxes  depicted  in  Figure  B 1  may  be  expressed  as 

$wa=  PwCCp  (rrrsw)  +  (rwL-rf)  lw  (bd 

hi 

in  which  the  second  term  on  the  right  side  refers  to  heat  conducted  through  the  body  of  lock 
water,  Tsw  is  the  surface  temperature  of  water  ( Tsw  >Tf)  and  Tw l  is  the  water  temperature  in 
the  lower  layer  of  the  lock  chamber.  As  an  approximation,  one  might  assume  Tw l  =  4°C. 
Using 


<t»wa  =  Wwa(TSw-ra)LW 


eq  B1  becomes 

q  _  Tj  *^W(rwL- 

PwCp(rg-rsw)  /t,pwCp(rg-rsw) 

_  LW  fj  (Tsw~  ^~a)  ky/  (^wL-  ^*f) 
Pw^p  (^g —  7"sw)  (T'g —  Tsw) 


(B2) 


(B3) 


It  is  assumed  in  the  foregoing  that  groundwater  discharge  cools  from  7g  to  Tjw.  the  surface 
water  temperature,  before  descending  to  the  lower  part  of  a  lock  chamber. 

In  eq  B3  the  thickness  and  temperature  of  the  lower  layer  are  not  constant.  Layer  thickness 
hi  increases  with  time: 


hx  = 


o_ 

LW 


ft 

Qdi 
.  o 


(B4) 


with  boundary  conditions:  ht  =  0  at  t  =  0;  and  h{-D  when  l  =  LWDIQ.  Layer  temperature 
decreases  at  a  rate 


<frw  _  1 
df  ht 


-^-(TwL-Twf) 

PwCp 


=  _2_ 

LW 


(Fg  —  Tsv/) . 


(B5) 
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